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ABSTRACT 


The blade-element performance of a small transonic 
compressor rotor operating at high subsonic speeds is 
reported. Measurements were made downstream using a small 
calibrated pneumatic probe, with simultaneous measurements 
of total pressure, velocity, turning angle, pitch angle, 
and temperature. In addition, a program is presented to 
calculate blade-element performance based on empirical data 


contained in NASA SP-36, 
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I. INTRODUCTION 


The study Periibed here 1s part of an ongoing program 
at the Naval Postgraduate School to determine the perform- 
ance of a transonic compressor stage. The overall goal of 
the program is to develop the knowledge necessary to design 
compressors with multiple transonic stages which are effi- 
Cient, light in weight, and free of flutter problems. It 
has long been appreciated that increasing the wheel speed 
and relative air velocities in axial compressors would re- 
sult in fewer stages for a given energy addition. However, 
in order to achieve high overall efficiencies in multiple 
axial stages, the ability to predict the behaviour of a 
Single stage correctly must be shown. The purpose of this 
study was to measure the blade-element performance of a 
transonic compressor rotor operating at high subsonic speeds 
and to compare the measurements with the performance calcu- 
lated according to the method given in Reference li. 

The transonic compressor and drive turbine were de- 
Signed by Dr. M. H. Vavra., The program is supported by the 
Naval Air System Command, Code 310, through the office of 
Dr. H. J. Muel ier, 

To determine the blade-element performance of a com- 
pressor blade row, it is necessary to evaluate radial veloce- 
ity and angle ence intemand ousevot the bikade row, 
In a transonic compressor, at or near design speed, the 


unsteady nature of the flows with large pressure amplitudes 
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complicates the problem. In addition, the possibility that 
wakes might induce vibrations in the highly stressed blades 
prohibit the use of probes upstream of the rotor. Measure- 
ments can be made downstream of the rotor; however, Since 
the flow area is small, boundary effects on probe measure- 
ments can be significant. 

In order to carry out the present study, a small com- 
bination temperature-pneumatic probe was designed by Dodge 
(Ref, 2) specifically to be used immediately downstream of 
the rotor. In previous work, the flow upstream of the rotor 
was determined at different flow rates (Ref. 3). It was 
left to the present work to represent the inlet flow field 
analytically, to calibrate the combination probe for bound- 
ary effects, and both to measure and to calculate the rotor 
blade element performance. 

In this report the transonic compressor and test rig 
are first described in Section II. The method of approach 
and program of measurements are described in Section III. 

The results of the measurements are then described and 
discussed ia weetion IV, and Conclusionsware given in Sec-~ 
Tron’, “Mnmonder to maintain a clear presentation of the 
relevant results, the details of the calculations and of 


the calibration measurements are given in separate Appendices. 
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II. THE TRANSONIC COMPRESSOR 


x. DESCRIFTION 

A drawing of the compressor test rig is shown in Fig. 1. 
The test rig consists of an air-turbine drive unit and an 
induction section containing a filter, a throttle, a set- 
tling chamber and a flow measuring nozzle. The air turbine 
is designed to deliver 450 horsepower at 30,000 RPM, The 
relative blade tip Mach number of the rotor at design point 
is 1.5. The compressor flow rate is controlled by an elec- 
tric-hydraulic rotating throttle plate. Power is provided 
by a compressed air supply from an Allis-Chalmers multi- 
stage axial compressor. A complete description of the 
Turbopropulsion Laboratory and test facilities is contained 
in Ref. 4% <A drawing of the transonic compressor is shown 
Dy I ta 2 - 

The blade design for the transonic compressor rotor was 
derived from circular-are elements wrapped on six conical 
surfaces. The outer surface (at the case wall) is a cylin- 
der. The blading geometry is given in Fig. 3. The blade 
pressure side is flat and the suction side is a circular- 
amc. “uae Soltvaicvy is fairly constant with a value of 1,35 
at the mid-streamline. The outer diameter of the rotor 
blade is 11 inches with a channel height of 1.92 inches 


downstream and a radius ratio at the inlet face of 0.5. 


16 





B. INSTRUMENTATION 

The transonic compressor iS equipped with stationary 
probes to provide routine performance data. Triple redun- 
dancy is provided in the measurement of the power; either 
the mass flow rate and temperature drop through the drive 
turbine, or the mass flow rate and temperature increase 
througn the compressor, or the speed and the measurement of 
torque can be used. The torque is measured by strain gauges 
on flexures against which the stator assembly is free to 
rotate. Flow straignteners before and after the stage en- 
Sure that the stator and rotor torques are equal. The 
method based on torque and speed is the most accurate since 
it is independent of flow profiles. Static pressure taps 
were installed on both outer and inner walls of the com- 
pressor casing. The rotative speed of the test rotor was 
determined by a magnetic flux cutter and electronic counter, 
The fixed instrumentation will not be described further, 
It was used in the present tests only to provide the refer- 
ence conditions for the survey measurements, The compressor 
cylindrical casing was modified to accommodate a variety 
of probes and sensors which may be required in the program 
of steady state and unsteady flow pressure measurements, 
The casing can ve rotated to position probes in the periph- 
emal dawection. 

The rotor exit survey measurements reported here were 
obtained using a combination probe (Fig. 4) designed by 


Lt. F. J. Dodge (Ref. 2). The combination probe, following 
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calibration, was used to determine total pressure rise, to- 
tal temperature rise, and flow velocity in magnitude and 
direction. In addition, static pressure was measured using 
the wedge static probe shown in Fig. 5, which was based on. 

a NASA design provided vy Mr. Lloyd N. Krause of NASA Lewis 
Research Center. Each probe had two “static pressure™ sen- 
sors whtch, if the probes were rotated until the indicated 
Static pressures were equal, allowed the probes to be aligned 
with the flow direction. It was determined that petter 
accuracy could be attained, however, by aligning both probes 
to the flow direction sensed by the combination probe. This 
was because the sensors on tne wedge probe were separated 

in the radial direction. A chromel-constantan fine-wire 
thermocouple in the combination probe was used in conjunction 
With a similar sensor mounted upstream of the rotor to 
measure, differentially, the total temperature rise, 

Radial surveys were taken at a fixed peripheral case 
angle py manually positioning the prove in the flow using a 
United Sensor traverse unit, manually rotating the combina- 
tion probe to balance indicated static pressures, and ad- 
justing the wedge probe to the same angle, All pressures 
were measured by a calibrated transducer ina Scanivalve 
arrangement. (A capability exists to remotely control probe 
movement and balancing through the use of a programmable 
microvrocessor and a probe traversing mechanism designed by 
Lt, D. D. Patton (Ref. 5). Remote operations will be needed 


as higher operational speeds are developed.) 
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C. DATA ACQUISITION AND REDUCTION 

Pressures, temperatures, rotational speed and other data 
were scanned and punched onto paper tape by a B & F Model 
24133 solid state data acquisition system. A complete de- 
scription of this system is contained in Ref. 4. The 
punched tape data was processed using Hewlett-Packard Model 
9830A programmable calculator data system. The data was 
transferred to a Model 9867B mass memory disk unit for stor- 
age prior to data reduction. The data facilities are shown 
in Mas. 5. 

A program, designated "WHO01" was written to transfer 
raw data from tapes to the mass memory. <A second program 
"HAWK1" was written to calculate temperature, corrected 
non-dimensional velocity, turning angle, total pressure, 
and static pressure. Two other programs, used in Sequence 
("REPROL™ and "RASC™) completed the data reduction process, 
Plotting routines developed by Dr. R. FP. Shreeve were 
adapted to plot distributions of Mach number, pressure rise, 
temperature rise, flow angle, and loss coefficients, using 
the X-Y plotter. 

In reducing the probe measurements to obtain velocities 
and loss data, the thermocouple calibration reported in 
Ref, 2 was used. However, new pneumatic calibrations were 
performed and corrections were derived as described in 


Appendix B, 
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iit, METHOD OF APPROACH 


A. METHOD 

The usual approach to the design of axial compressors is 
described in Ref. 1. The method involves combining an axi- 
symmetric through-flow analysis, which prescribes the stream 
surfaces, with a two-dimensional cascade performance calcu- 
lation which determines the required blading geometry on 
those surfaces. The blading is generally layed out on coni- 
cal approximations of the prescribed stream surfaces, 

When a new machine is built, the opportunity is there 
to compare the actual blade element performance established 
by measurement with that which can be calculated by applying 
the design approach. The comparison can be used to derive 
corrections for use in future designs, and it is particularly 
useful to obtain data for new blade profiles. 

The aaeeeae compressor investigated here was designed 
using the above approach, and the resulting blading is shown 
ieee, 3. Lt as to be noted that the rotor blading is not 
typical in design. In order to establish the blade element 
performance of a rotor experimentally, measurements should 
be made of the flow into and ait of the rotor on the stream 
Surfaces on which the blading profiles are prescribed, How- 
ever, it is desirable not to operate at high wheel speeds 
with strong wakes from upstream probes, and so in the pre- 
sent work a method of prescribing the inlet flow field from 
measurements of only total flow rate and wheel speed was 


derived. The method is given in the following section, 
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The flow downstream of the rotor was determined using two 
calibrated probes, The calibration and corrections re- 
quired for the probes are described in Appendix B, From 
the rotor exit flow measurements, the blade element loss 
coefficient and the local deviation angle were calculated. 
Finally, as described in Appendix A. the method of 
calculating deviation angle and loss coefficient given in 
Chapter 7 of Reference 1 was programmed for the Hewlett- 
Packard Model 9830A calculator so that these quantities 
could be computed for four blade element surfaces at any 


given compressor flow rate and operating speed, 


B. INLET PROW FIELD REPRESENTATION 

The flow field ahead of the rotor was measured in early 
tests using a calibrated pneumatic probe (Ref. 3). The 
measurements were compared with the results of a compress- 
ible finite difference calculation reported in Reference 6, 
and blockage factors were derived to relate the inviscid 
calculations to the measured flow field. 

The method adopted nere was based only on Anderson's 
measured data, As described in Appendix C a polynomial re- 
lation was found to correlate the measured velocity distri- 
butions away from the hub and tip regions, The flow function 
in this region was therefore calculated using Eq. C(3) and 


iw C(7). 
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C. Us® OF FLOW PROBES 

Pressure probes are designed to exploit the distribution 
of pressure which occurs over a body immersed in a moving 
fluid. These pressure variations depend mainly on the flow 
velocity so that with a carefully selected choice of body 


shape, having multiple sensors,a probe can be calibrated to 


) 
determine flow velocity. The relationship between pressure 
and velocity must be established over a range of Mach nun- 
bers. In addition to flow velocity in magnitude and direc- 
tion, the distributions of stagnation temperature, and total 
and static pressure are usually required. The probe designed 
by Dodge and reported in Ref, 2 is a small lightweight probe 
of aerodynamic design which allows these quantities to be 
determined at any position in the flow, 

In order to provide an independent measurement to com- 
pare with the combination probe, the static pressure distri- 
bution was measured with a second probe, and static pressures 
were also recorded at hub and tip wall taps. The wedge 
Static probe shown in Fig. 5 has the favorable characteristic 
that when it is aligned with the flow in the yaw plane (by 
balancing the pressures measured on the wedge surfaces) the 
indicated static pressure is relatively insensitive to angle 
variation in the pitch plane. As reported in Ref. 7 this 
type of probe has been used successfully up to a Mach number 
of 1.6. Further, for measurements in the transonic regime 
the probe has favorable characteristics. If the leading 


edge is sharp and the wedge is small, Ref. 7 states that the 
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calibration is also independent of Mach number. This was 


found to be the case for the probe used in the present work. 


D. PROGRAM OF MEASUREMENTS 

The combination probe and the wedge static probe were 
mounted at the exit of the rotor with a peripheral separa- 
tion of 120 degrees, Four Surveys were conducted at 50 
percent of design speed, at throttle conditions varying 
from near surge to maximum flow rate. Case angle for the 
first survey was 345 degrees. Subsequent surveys were con- 
ducted at a case angle of 357 degrees because it was deter- 
mined that peripheral variations induced by the stator 
blading were felt least at this setting. RPM was the same 
(e5,2902 390) for @ll four runs, 

The procedure in each survey was to move the probes in- 
ward Simultaneously, balance the combination probe to sense 
the flow angle, and then set the wedge probe to the measured 
angle, A data scan was taken, and the procedure was repeated 
at 0.1 inch intervals. Because of the hooked geometry of 
the wedge probe, measurements with this probe of static 
pressure could not be made closer than .6 inches to the 
outer wall boundary. Static pressures were derived for the 
complete passage through the use of the combination probe 
calibration, Wall static pressure measurements at the nub 
and tip were also recorded, Yaw angle and linear displace- 
ment were recorded manually. All other data were automati- 
cally recorded by the data recording system. Pressures were 


recorded to an accuracy of .1 inches of water. 
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V. RESULTS AND DISCUSSION 


A, TIME-AVERAGED FLOW FIELD 

The distributions of total pressure rise, total temper- 
ature rise, Mach number and flow angle measured using the 
combination probe downstream of the rotor are shown in 
Pras. 7, 6;°9, amd To. 

The total pressure rise was nearly constant with radius, 
whereas tne temperature rise was indicated to be larger 
toward the tip. The larger indicated temperature rise near 
the case wall might be caused by higher temperatures in the 
incoming case wall boundary layer. Both the Mach number 
and the flow angle increase toward the hub, so tnat the 
througn-flow velocity tends to remain constant. The results 
for the flow angle are qualitatively consistent; however, a 
large decrease can be observed in the magnitude of the flow 
angle for the open throttle case (Fig. 8). No definitive 
explanation was found, but as will be shown later this 
measurement leads to inconsistencies in the correlation of 
olade element performance. 

The velocity and pitch angle distributions, which were 
derived from the probe calibration and corrected for bound- 
any trees (Appendix B), are shown in Fig. 11 and Fig. 12. 
The pitcn angle shows a smooth variation between the values 
of the wall slopes at the hub and tip boundaries. It is 
Shown in Appendix B that the pitch angle is determined to 
only accuracy, So that the variations in pitch angie 


indicated in Fig. 12 might be exaggerated. 
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A Sensitive test of the accuracy of the combination probe 
measurements- was to examine the static pressure distribution 
inferred from the measured velocity and probe Ne pres- 
Sure measurement, Fig. 13 shows the distribution of static 
pressure calculated in this wax, Also included in Fig. 13 
are the values of static pressure measured at the hub and 
tip wall static taps. It can be seen that agreement near 
the case wall was excellent. The comparison near the hub 
wall was not as zood; however, very Similar results were 
obtained from the wedge probe as shown in Fig. 14. While 
no readings were possible near the outer wall using the 
wedge probe, the two probes gave very Similar results near 
the hub wall, and gave almost identical results for the sur- 
vey at open throttle. 

As shown in Appendix B, the effect of flow boundaries on 
the readings of the combination probe is small. The magni- 
tude of the correction which was included in the above re- 
Gules Can be seen in Fig, 15. The correction is less than 
1% in velocity over the useful range of the survey. 

Velocity diagrams near the mean radius (0.55 of the 


channel height) are shown in the four parts of Fig. 16. 


B, DEVIATION ANGLES 

The design procedure for calculating the deviation 
angles for a given blade profile, with corrections for 
three-dimensional effects, is described in Appendix A, 


Blade profiles are known on four conical surfaces in the 
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flow field, on which the blading dimensions are those given 
on Fie. 3. Using relative flow angle Bo from the velocity 

diagram data (Fig. 16 for example), the flow deviation angle 
($8) can be calculated from the stagger angle (¥) and camber 


/ 
angle (9) using the relation 


& = @,-%+4 


/ 


The results for the four probe surveys are shown in 
Fig, 17.-1n comparison with the results of calculations 
Following Ref, 1. 

Departures from values predicted using Ref, 1 are quite 
large. However, the design procedure strictly applies only 
to the case where the operation is at the point of minimum 
loss for the given blade element. <A closer correlation 
might be expected as tne higher operating levels are achieved 
and closer to design conditions are attained. (Ref. 1 states 
that at transonic inlet Mach number levels, the point of 
minimum loss should constitute the desired design setting. ) 
It is noted that the through-flow velocity is considerably 
increased through the rotor at the 50% design speed setting, 
and this would effect the deviation angle. No explanation 
was found for the indication of negative values of deviation 
angle near the hub, although a careful check was made on the 
accuracy of the probe yaw angle measurements and on the ac- 
curacy of the data reduction program. It must tentatively 
be concluded that there are strong three-dimensional effects 
present which can lead to the indication (locally) of nega- 
tive deviation angles. More data must be obtained to investi- 


gate tnis speculation. 
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C,. LOSSES 

Loss coefficients for the flow over the blade elements, 
assuming that the stream surfaces were the conical surfaces 
given in Fig. 3, were calculated from the probe measurements 
uSing the methods described in Appendix D. The calculated 
loss coefficient (W) was consistent with the definition 
of loss coefficient used in Ref. 1. Of the two methods 
described in Appendix D, Method 1 uses both the stagnation 
pressure rise and stagnation temperature rise measurements 
of the probe. Method 2 does not use the stagnation temper- 
ature rise measurement, but uses a temperature rise calcu- 
lated from the velocity and angle measurements, through 
Euler's turbine equation. 

The results are snown in Fig. 18. Good agreement was 
obtained between the two methods of calculating the loss 
coefficient, however negative values were calculated for the 
inner stream surface, It should be noted that all measured 
data and all blade properties varied smoothly with radius, 
so that smooth curves between the individual data values 
could be drawn, The distributions are in qualitative agree- 
ment with earlier measurements, and it is noted that nega- 
tive values of losses were previously reported by Dodge 
(Ref. 2) from measurements made with the same combination 
probe, 

A negative value of the loss coefficient can result only 
from a measurement of stagnation pressure rise which is too 


high, or a measurement of stagnation temperature rise which 
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is too low. This can be seen if Eq. D(13) is rewritten in 


the form 


[1 rn xo] Y 4-1) to 


die 
ye 


where Wis the ratio of total temperatures, (Ty,/T+,). 


| E - 4] 


[1 : 2 a 


Clearly, since Xy5 and Xy, are positive,(o is negative only 
if yet > (Pt+o/P+,)- pinee the prediction of Method 1 and 
Method 2 agree, and because a favorable comparison with ear- 
lier unreported measurements of temperature rise was made, 
it is tentatively concluded that the impact pressure rise 
measurement is too large. <A positive error in impact pres- 
sure could result from the pneumatic integration of the peri- 
odic unsteadiness of the flow from the rotor. All parts of 
the probe are the open ends of continuous tubing. It is 
known that pneumatic probes of different designs integrate 
Geencwently in periodic flow conditions. An investigation 
of the pneumatic averaging characteristics of the probe is 
required in order to correct the measurements of losses, 

It should be noted that the velocity and pitch angle calcu- 
lated from the time-averaged probe readings should not be 

as sensitive to the averaging behaviour of the probe. This 
1s because all parts are made from similar tubing and will 
average ina Similar way. The loss coefficient, however, 

is very sensitive to errors in impact pressure measurement 


alone, 
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While the magnitude of the loss coefficients is clearly 
in ercemygei teks Of intemest to attempt to commelate the 
measurements on the basis of diffusion factor as defined in 
foie lw Theseresult of this attempt is shown in Fig. 19. The 


Gaeewenem acter (D) was calculated from the definition 


Wo RWy, + RoWy 
aoe. 4 2" u2 


“LO (R, + RS) W, 
where W denotes the relative velocity, R is the radius, O is 
the solidity and subscripts 1 and 2 denote stations upstream 
and downstream of the rotor, re is the peripheral component 
of the relative velocity. The “total loss parameter", 
(cos 8 ./20 ) is consistent with that shown in Fig. 203 
of Ref. 1 and was calculated from probe measurements, 

It can be seen in Fig. 19 that a reasonable correlation 
of the loss coefficient data on the basis of diffusion factor 
was obtained for three of the four sets of data. The data 
for full open throttle departed significantly from the data 
irom three other test conditions. It should now be recalled 
that the flow angles measured for this test were also incon- 
Sistent. More data is again needed, and a correction must 
be devised for unsteady flow effects on the probe measure- 


ments before a complete correlation can be obtained, 
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VI. CONCLUSIONS 


The attempt was made here to determine the blade element 
performance of an "impulse" rotor designed for transonic flow 
(design point) but operating at high subsonic (off-design) 
conditions. The results were preliminary, but served to 
define the problems involved in obtaining necessary design 
data from empirical correlations of measurements made in 
actual machines. 

The results obtained depended on the performance of a 
particular combination probe when ret downstream of the 
rotor. The following conclusions were drawn concerning the 
probe measurements : 

(i) The effect of the presence of flow boundaries on the 
probe measurements was small, but corrections were derived 
and apolied in the present work. Static vressure derived 
from the probe measurements, static pressure from wedge probe 
measurements and wall static pressures were in reasonable 
agreement, 

(ii) The temperature rise across the rotor measured by 
the probe, using the vrobde calibration, was consistent with 
the rise deduced from Euler's turbine equation, All indica- 
tions were that the probe provided an accurate measurement 
of the temperature rise. 

(iii) The impact pressure measured by the probe in the 
unsteady flow downstream of the rotor was probably too high. 


An investigation of the unsteady response of the tubing used 
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in the probe design is required in order to correct the 
data obtained so far. 

(iv) The flow velocity derived from the probe measure- 
ments was probably not as sensitive to unsteady flow effects 
Since all the pressure sensors on the probe were constructed 
from the same tubing. The derived velocity depends mainly 
on differences between the pressures registered by the sen- 
sors. In steady flow the velocity was determined to better 
than 1.0% accuracy, to nearly sonic speeds. 

(v) The pitch angle was determined by the same probe to 
within 13°, 

(vi) The wedge probe provided a good measurement of 
Static pressure in the compressor when rotated to a pre- 
measured flow angle. 

As a result of (iii), the loss coefficients determined 
from probe measurements were too small, becoming negative 
Omer tie Diadine near to the hub. A correction to the im- 
pact pressure measurements and more data are needed before 
a useful correlation of the losses can be obtained. However, 
1t was found that the present results from three tests and 
four different blade-element surfaces correlated reasonably 
well in terms of the NASA diffusion factor. 

The deviation angles measured in the present tests de- 
parted significantly from the values predicted using NASA 
SP-36,. No explanation was found for the negative values 
measured near to the hub, although the possibility exists 


that they could be caused by three-dimensional effects in the 
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flow through the rotor. More data are needed to confirm the 
trends in the deviation angles, Agreement with the predicted 
values was not expected since the blade profiles were unusual, 


and were operating well away from the design condition. 
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APPENDIX As BLADE ELEMENT PERFORMANCE PREDICTION 





Al. INTRODUCTION 

Flow through the blading of an axial-flow compressor is 
an extremely complicated three-dimensional process, As a 
result, Simplified design models have been devised that per- 
mit the use of two-dimensional techniques. The approach is 
described fully in Reference 1. Involved is first the cal- 
culation of the “through-flow” assuming axial symmetry, then 
the calculation of the “blade-element flow" assuming approx- 
imately two-dimensional behaviour. The latter calculation, 
based on Chapters 6 and 7 of Reference 1, is the subject of 
mae SeCHion. 

Five quantities shown in Figure Al determine the aero- 
dynamic behaviour of a cascade in ideal flow. Three of these 
define the geometry of the blade; they are, the shape of the 
blade, the ondlentation of the blades with respect to the cas- 
cade axis, and the solidity. The fourth quantity identifies 
the direction of flow ahead of the cascade and finally, the 
Mach number must be known. When these quantities are speci- 
fied, the inviscid flow may be determined. 

Real cascade flows however, experience many other in- 
fluences including friction, turbulence, three-dimensional 
and unsteady effects, and boundary layer interaction. The 
procedure outlined in Reference i correlates cascade and com- 
pressor blade element performance data. Various correction 
factors are introduced which account empirically for the 


effects which can not be predicted analytically. The result 
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msi a method of calculating the ee ioemeinGs of compressor 
blading. The method can be applied to both rotors and stators 
over a fairly wide range of Mach numbers. The method is 
limited however to the calculation of performance close to 
design point, 

NASA SP-36 has become a reference text in the axial 
compressor industry. A brief summary of the method is pre- 
sented here, followed by a computerized solution which uses 
the extensive correlations of data presented in this manual, 

The main thrust was to develop a program to calculate 
the blade element performance of a blade design of arbitrary 
geometry. The program would provide a very powerful design 
tool, and would provide the basis for a program to predict 
the performance of a transonic compressor rotor. Reference 6 
reports an extension of this was used with considerable 


Success up to a Mach number of about 1.5. 


A2. METHOD OF BLADE ELEMENT PREDICTION 

A2,1 Two-Dimensional Incidence Angle. Any cascade of 
aerodynamic blading behaves qualitatively as shown in Fig- 
ure A2 where & i is the “range” of performance and the mid- 
point of the curve is defined as the reference minimum loss 


incidence angle (i ). For blading with a given camber 


wer. 
(S'), thickness. distribution, and maximum thickness ratio 
(t/c), the solidity (©) and blade angle (4%) complete the 
definition of the cascade geometry. Given@ , and the inlet 
velocity the method in Reference 1 allows the determination 


Minimum loss incidence angle, the deviation angle, and the 


loss coefficient. 


Dil 


The reference incidence angle for any blade element can 
be calculated from correlations of two-dimensional cascade 
measurements, Two-dimensional reference incidence angle for 
10% thick blades is given as 

is_p = Ks(igdio + ng' A(1) 
where K, is a function of blade thickness distribution and 
(igndso is the zero-camber incidence angle (function of G4 
and ), nis a slope factor given by 

ee ie i ie/ id oup A(2) 
which is also a function of @, and Cue. 

in Ref. 1, (indi iis “Civen@ins Fists 7 and nis "einen 
tn Figs. 1338, Ks has the value 0.7 for double circular-are 
aerofoil shaves, 

A2.2 Two-Dimensional Deviation Angle. As in the analy- 
Sis of reference incidence angle, the correlation of blade- 
element deviation angle is expressed in terms of blade 


geometry by 


Seep = Xe ($010 +{m /o "|p A(3) 
where Keis a function of blade type (0.7 for double circular 
arc). 4 ele is the zero-camber deviation angle (function 
of 8, andO), m is a function of @,, for the different basic 
Seer distribution, and b iS an exponent that is alsoa 
iwmetion of e 1° 

Since deviation angle varies somewhat with incidence 


angle, a corrected deviation angle is defined as 
O 
Oo _ O m es dg 
Yen : “ (30 10 + 26 P FT, - Aju) ee op 
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where (a$/d.)5_p is the slope of the two-dimensional varia- 
tion of deviation angle with incidence angle at reference 
incidence, and Le is the actual incidence angle. 
Camber and turning angle are related by the equation: 
p= 61 -R2tari 
or 
AQ= gr +i -§ A(5) 
Once the two-dimensional deviation and incidence angles have 
been established, the turning angle can be computed. 
A2.3 Corrections for Three-Dimensional Effects, The next 
step in the procedure is to correct reference deviation and 
reference incidence angles for three-dimensional effects, 
This is accomplished by defining the three-dimensional refer- 
ence angles 


ie? 5b = leap tte gapdee BN) 


and 


(Sref3-p ~$2-p + (8 - $2-p) A(7) 


where (i - i,_p) and (§ -§5_p) are functions of Mach number. 
In summary, the desired turning angle is given by Eq. 


260 “of Ref, 1. 


4@ =p [1 -3+n]- (io - See) 2g ~ 2305) 


ag° 
[1.0 - (F) op | - Kl 0a Ke ($0) 10 A(8) 


A2.4 Loss Coefficient, The loss coefficient for a blade 


element can be calculated from the correlation given in 


De 





Fig. 203 of Ref. 1. Here the “total loss parameter" is cor- 
related in terms of "Diffusion Factor", The calculation of 


losses is given in Appendix D. 


A3. POLYNOMIAL REPRESENTATION OF CORRELATIONS 

All values on the right side of Eq. A(8) may be deter- 
mined using data given in Figures in Ref, 1 once the geometry 
of the blade has been given and the inlet velocity calculated, 

Fach of the unknown values of Eq. A(8) are functions of 
blading geometry and air inlet angle or blading geometry and 
Mach number, and are described by well-behaved families of 
etevee, bie, 166 of Ref. 1 has been chosen to illustrate 
the process of computerizing the solution, which involved re- 
ducing the data correlations given in figures to analytical 
expressions, 

Eleven values of air inlet angle (6,) were chosen for 
eeen solidity oe in Fig. 168. The values were plotted 
as slope factor (m) vs. air inlet angle (8) for constant 
solidity. The data points for each solidity were fit by a 
third order polynomial curve using the Hewlett-Packard 9830A 
calculator and X-Y plotter. For simplicity, a fifth order 
plotting routine was always used with the higher order co- 
efficients equal to zero if a lower order resulted in a more 
aperopriase fit. The plot and polynomial curve fits for 
each solidity are shown in Fig. A3. The polynomial curve 
fit for each set of data points gave six soefficients for 
egen clive. Thus for each m vs. Se, plot, the six coeffi- 


cients Ags Ayneseeerde were calculated for the polynomial 
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expression for m as a function ener 


(a 


m= A eye, * 2 


: a” + ° in A, 8," + A 6,7 A(11) 

at each value of solidity with eC, converted to radians to 

minimize the error in higher order terms, The coefficients 

Ao» Ayreveeesde derived from the curve fits of Fig. A3 were 

plotted against each of twelve values of solidity as follows: 
Ay (I) vs. o (I)........12points 

Ween SO XT). ccs eee e Lepoints 

Y 2) 
hs Or) weme ae (1).5...-..12points 


This produced curves of the form 


— 
il 


3 yo. 
0 Bao + ao + ano a Azo tayo + aug? 


i! 1 ie? al 
‘me cre + alo + 856 + ayo? 


A 1 0 
A(13) 
A 


he eee 


: a? ee aio + a, 2 oF a ed a arg" 4 a 
which resulted in polynomial curve fits shown in Fig. A4. 
In this case the higher order coefficients An and Ag are set 
equal to zero. When all the coefficients in Eq. A(13) are 
known, the slope factor m, for any given oO and B45 can be 
calculated from: 

m = (a, yee ls ano” + 2x67 + ano + ant A(14) 
Lene es 


1 4 
+ (ao +a 16 +46 + a5 6 + a,°6 + a. = Ja 


a 


Sit(ay? + ay + os ys 767 + ayo" + 2s "5784 ; 
This procedure resulted in a Single equation (Eq. A(14))for 


m in which the coefficients were analytic functions of 6,. 


61 


’? 





Thescoetficients are given in Tables Al and A2. Due to the 
well-behaved nature of the curves, the slope factor m could 
be @eiculated to well wathin the accuracy of the omen 
data contained in Ref. 1, 

All data necessary to predict performance for circular- 
arc blading was reduced to polynomials using a similar 
procedure, The polynomials were included in a computer 
program to calculate blade element performance for specified 


ifeut geometry and flow conditions, 


A4, COMPUTER PROGRAM 

A4.1 Description. Program "NASA36" was written for the 
Hewlett-Packard 9830A calculator to calculate the incidence 
and deviation angles for a rotor with meridional inlet ve- 
locity. The program follows the method of Ref, 1 through 
the equations given in Section AZ. Where a value from a 
hasurme in Ref, 1 is required, it is calculated using a poly- 
nomial approximation to the data as described in Section A3. 

The program is self-explanatory through the listing 

given in Table A3. The symbols are defined in Table AH, 
The program performs the blade element calculations for a 
Specified compressor operating condition, The inputs are 
referred flow rate and speed, and temperatures upstream and 
downstream of the rotor. Through the method given in Appen- 
dix C, these inputs define the inlet velocity as a function 
or the radial Seuation. It should be noted therefore that 


the present version of the program includes a description 
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of the transonic compressor inlet flow-field, while it 
allows changes in blading geometry to be input. 

A4,2 OQOveration. Detailed information on the operating 
procedures for the HP-9830 calculator, mass memory and as- 
Sociated equipment are found in Ref. 9 and Ref, 10. The 
following is a step-by-step operating procedure for the 
pmecmam called “NASA36". 

1, After loading "NASA36" into the Hewlett-Packard 
9830A calculator press “Run" and "Execute", 

2. The calculator will then display a series of inputs 
required to determine the geometry of the blading. These 
are shown in Fig. 3 and are: 

a) Stagger angle 

b) Chord with sharp leading and trailing edges 
c) Tip radius 

d)-Hub radius 

e) Diameter of cone through stacking radius 

f) Maximum Profile Thickness 

eo) Radius of Profile Contour 

h) Blade camber angle 

i) Profile chord 

j) Angle of cone thru stacking radius 

3. The operating parameters of the compressor rotor 
are now required. These are: 

a) Revolutions per minute 
b) Total temperature upstream of rotor 


c) Temperature downstream of rotor 
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ad) Referred flow rate 
Outputs include deviation and incidence angles, Mach number 


before and after the rotor, and referred flow function (@). 
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APPENDIX B: CALIBRATION AND CORRECTION OF PNEUMATIC PROBES 


Bi, FREE-JET APPARATUS 
The free-jet apparatus used to calibrate probes is 
shown in Fig. Bl. The device consisted of a 8 inch diameter 
Pipe with the air supplied by the laboratory's central 
Allis-Chalmers Compressor, A 4.2 inch converging nozzle 
was used to increase velocity and attain Mach numbers up 
to about .9. To determine if the exiting velocity was uni- 
form and undisturbed by the probe mounting arrangement, 
Surveys were conducted with a Prandtl probe over a range 
of 5 Mach numbers at two axial locations. osSurveys were 
also conducted with the same probe mounted and traversed 
from opposite sides of the jet. The results shown in 
Hie. Be andde@ated that the velocmey profiles were uniform 
over the core of the jet extending to a diameter of 3.8 
inches. As suggested in Ref. 11, the non-dimensional velo- 
city (X) appearing in Fig. B2 is defined oy 
TIONS Ie a) 
a! sosiade: O : 
Where V, = 2c 7, = 10%.62/T, when T, is in R. Following 


Ref, 11 further, the relationships for stagnation pressure, 


total temperature, and Mach number are, in terms of X given 


by 
pT — : 
= 1 - xX B(2) 
Poe. - x2) M01) B(3) 
= 
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VYRT aa 1 - X il 

The values of X in Fig. B2 were calculated using the meas- 
urement of stagnation pressure and the assumption of atmos- 


pheric static pressure in Eq. B(3). 


B2. CALIBRATION OF THE COMBINATION PROBE 

Calibration of the pneumatic probe reported in Ref, 2 
can be represented by defining two functions of the non- 
dimensional velocity, X, and the pitch angle, g. If P, 
and Po /3 are indicated stagnation and indicated static pres- 


Sure of tne probe respectively, then define 


ls Mali 
Sam -— 4, B(5) 


Fe = i 
P, = P 
By = Fp (x) = + aie 


to be calibration functions when the probe is rotated to 
balance P., = Pa. Pi is the fourth, off axis pressure meas- 
ured by the vrobe. 

The complete calibration procedure is given in detail 
in Ref. 3. The prove was recalibrated to confirm the re- 
Sults given in Ref. 2 before investigating corrections for 
boundary effects. A range of Mach numbers (X = .13 to .4) 
and vitch angles (-20° to 20°) were used, With true static 
pressure taken to be atmospheric and true impact pressure 


that measured by a reference Prandtl probe, the value of 


A @e@s calculated from Eq. B(3). As in Ref, 2, polynomial 


83 


- —_- eel 








= 
ke 
ie ee <—_> _ en 
_ 
- = ~~ 
> 
Sr a = 
= _ —— 





representations were found for the calibration equations 


written in the form 


X= X(B, @) B(7) 
and 

Fp = Fp (X,P) B(8) 
Where F,., and Fy are the ratios of probe pressure measure- 
ments defined in Eq. B(5) and Eq. B(6). it is moved char 
Eq. B(6) departs from the definition used by Dodge in Ref, 
2. The present method was found to give more accurate re- 
sults. The coefficients in Eq. B(7) and Eq. B(8) derived 
from the calibration test results are given in Table Bi. 

As a test of the analytical representation of the 
calibration, the calibration data were used to calculate 
values of X and ®@to compare with the actual values. Table 
Be shows the results of this comparison for the velocity. 
Te error in X is generally less than 1% and only reaches 
nigher values at high positive or negative pitch angles 
( 16° or -8°) or at higher Mach numbers (.88). At this 
point the error increases to about 2%, 

Pie, B3 shows the error in the piten angle in tne 
mneamee of piten angle from ~4° +9 +12°, The error is seen 
to be 714° and is not a simple function of either pitch 
angle or Mach number (or X). ‘There appears to be a limit 
to the capability of resolving D. The error was larger at 


high negative and positive values of pitch angles but these 


were outside the intended range of uSe. 
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B83. BOUNDARY EFFECT ON COMBINATION PROBE MEASUREMENTS 

In application, the probe is not in plane flow away 
from all boundaries. The indicated velocity and pitch angle 
will depend upon several factors including boundary layer 
and stem immersion effects, The method of approach assumed 
here was that the probe always obeys the calibrations es- 
tablished and that errors are induced due to the perturbation 
of the probe itself on the flow near the boundary. Thus, 
in any flow the probe measures a velocity, X, and a pitch 
angle, (oF which are the sum of the actual and an “induced” 
component, 

Three surveys were conducted ina 10 inch pipe up to 
She maxamtm attainable Mach number of .55. The results of 
calculating velocity applying the free jet calibration given 
above are shown in Fig. B4. 

Introducing subscripts a(actual), i(induced), and 


m{measured), the velocity, X, may be represented by 


2 
—_ = i = 
as, 10k, = %, [p+ eh] X, (1 +€,) B(9) 


where the fractional ‘error' in X,€,, which is induced by 
the presence of the wall, is given by 


XX =- X 


~~ nh a 
Eee ——— B(10) 


He js determined from the probe calibration in the usual 
way. 4%, 1s the required, correct velocity, and€, is the 


error to be established by calibration tests, 
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The measured velocity is corrected using 


Xx 
= mm 
Ke el me Bigiaiy) 
x 


when the error€, is known. 
In the same manner actual pitch angle may be repre- 


Sented by 


P,.2=9-& B(12) 


In the 10 inch pipe tests, the actual velocity was 
calculated from the probe impact pressure and pipe wall 
Static pressure using Eq. B(3) and the actual pitch angle 
wasezero. Thus the “induced” error in X was cahculated 
uSing Eq. B(10) and the "induced" error in g was calculated 
same Bq. B(12). The two errors, €, and ey , were then ex- 
amined as a function of distance from the pipe walls ex- 
pressed as a number of ‘probe diameters; (the characteristic 
dimension of the probe was taken to be 0.1 =g/inches). It 
was found that the errors became negligible at about 9.5 
probe diameters from the boundary. In order to derive an 
expression of a suitable form, the errors were plotted as 
a function of (9.5-D), where D is the displacement as a 
number of probe diameters, The velocity error was then 
approximated by a polynomial of the form 
2 


é.. ~ A a: A, (9.5 rx D) ah A, (9.5 a D Bigis:) 


O 
where D is measured from the wall through which the probe 
was inserted. The result is shown in Fig. B5. The error 


was less than 2% at distances greater the 0.25 inches from 
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the wall, In order to extrapolate the correction to larger 

values of velocity, the error was represented as a function 
of X, The coefficients in Eq. B(13) were plotted as a func- 
tion of X, and the result is shown in Fig. Bé6, Ao» A, and. 

A, were seen to be almost linear functions of X and were 


therefore approximated by 


— 2 
Ay = By ~ Bi os > Box os 
Mm =c +¢6.x%..+c.x .” 
| 0 17.95 2° .95 
iD te Diem: Dex. ° 
2 O 17.95 lemme 


where X gs was the value of X at D= 9.5. Values of 
X(corrected) minus X(uncorrected) verses immersion depth are 
Shown in Fig. 37. Because of the aerodynamic design of the 
probe, the corrections necessary for boundary effects are 
small. 

Results for the induced pitch angle error are shown in 
Fig. B8. Since these errors were within the resolution 


capability of the probe (14°), no correction was made, 


B4, COMBINATION PROBE MEASUREMENTS IN STEADY SWIRLING FLOW 
Additional errors may be present when a pneumatic probe 
is immersed in a swirling flow. Fig. B9 shows the apparatus 
which was constructed to investigate this possibility. The 
apparatus consisted of an annulus with a channel height of 
14 inches (10 inches in outer diameter). A total of 40 
static wall readings were available and a provision was 


made for probe surveys at different axial stations. The 
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Swirl velocity was generated by a single stage of stator 
blading at the entrance to the annulus, Initial measure- 
ments suggested the possibility of boundary layer separation 
amd reverse flow near the hub. In an effort to correct this 
problem, a plate was constructed to decrease the channel 
height at the exit, isolating the pressure at the hub wall 
from the surrounding atmosphere, and olocking backflow, 
Results are Ziven in Fig. BlOwhich shows the static 
pressure derived from tne probe calibration. Agreement 
is very close near the hub with a moderate departure near 


tre tip at the highest Mach number, 


Bom CrbLemATION OF THE WEDGE PROBE 

A type of probe which is reported to give a good indi- 
cation of static vressure over a wide range of Mach number 
is the wedge static prove illustrated in Fig. 5. The probe 
nas a Meciet-shaved head with small pressure taps in eacn 
face of the wedge. In principle, if the leading edge is 
sharp and the wedge angle is small, the flow past the 
sensing orifice will be almost unperturbed, and the meas- 
ured pressure will be close to the static pressure in the 
undisturbed flow. 

The static pressure may thus be determined by 


-~ P, ) +P B(11) 


273 A 


where it is assumed the probe is aligned in the flow with 


2 Gs 


S A 


P, and P. equal, and P, is the reference atmospheric 


2 2 


oressure, 
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The wedge probe shown in Fig. 5 was traversed in both 
the free jet and the swirling annular flows through the 
Same range of Mach numbers reported for tne combination 
prove. However, the characteristics of the probe in 
swirling flow were not fully determined. It was found in 
analyzing the results of the tests in the swirl annulus 
that the radial gradient in the flow angle gave rise to 
probe readings which were in error when the probe was ro- 
tated to balance. This was explained by the difference in 
the radial position of the sensors on the two faces of the 
wedge, It was therefore necessary to use the combination 
probe to sense flow direction in the compressor measure- 


ments and to set the wedge to the measured angle, 
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APPENDIX C: INLET FLOW FIELD DETERMINATION 





Ci. METHOD OF APPROACH 

An important consideration in a transonic compressor 
installation is the velocity profile at the face of the 
rotor. The profile depends on the inlet shape and the mass 
flow rate. The following analysis is based on the assump- 
moon that a Close approximation to the true velocity may 
be calculated if the referred flow rate is known, 

Seven surveys were conducted by Anderson (Ref. 3) in 
whicn accurate calculation of the inlet velocity profiles 
was accomplished using a calibrated United Sensor Probe 
over a range of flow rates tabulated in Table Cl. 

Reference ii defines a “total flow function,” , given 


by 





cy = Rv = whe = F(X) = X(1 - x21 /(8-1) eC) 


The flow function provides the means of replacing density 
in the equation with the non-dimensional velocity (X) for 
internal compressible flows. For each survey conducted in 
wet. 3, © was calculated at each point and plotted versus 
the normalized displacement (hub to tip). The results are 
shown in Fig. Ci. Each curve was then fitted with a suit- 


able polynomial of the form 


a z 3 
o =a, + a ,H +t A oH a agit nae eee Ce) 


where H = normalized displacement (hub to tip) and Ass Anree 
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are determined by a polynomial regression curve fitting 
routine. The value of the total flow function at the cen- 
terline (Soy) was evaluated for each case. Values of 
255 for each survey were then calculated and plotted 
against H. The resulting curves shown in Fig. C2 were then 


approximated by a family of polynomials having the form 
oO = wo = Ao a AH a Co Ci) 


f 
which is then called the ‘normalized flow function, 


The referred flow rate is defined as 


pe. Se o (4) 
S 
mnere Tt, 
a= C(5) 
pet 
d 
an Pe, 
ae oe o(6) 
Ger 


and W = mass flow rate. 
% 
A’ plot of Boy versus W shown in Fig. C3 yielded the 


following linear relationship: 


Eo = dy + d,W C(7) 


Thus, Eq. C(3) and Eq. C(7) provide an analytical rep- 
resentation of the inlet flow profile in terms of the 


referred flow rate only. 
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Cee RPPELCATION 
Mee Gach run of the transonic compressor,W is calcu- 
lated and $ . is then evaluated from Eq. (7). ® can then 
be evaluated from Eq. C(3). 
- -X is then evaluated as the subsonic solution of the 


flow function given by Eq. C(1) using Newton's method as 


follows: 

1) Guess X = Ko 

2) Calculate 
Gort x2 yl a) c(8) 
dé _ 1 2 X 
dx ~ 22 + ~Wat) tales) | a 
x = & ee ¢(10) 

ax 


3) Now test for convergence of @, 
bet \?., - 8 , |<.00001 then the calculation is complete. 
If this requirement is not satisfied then the new 2, in 
Eq. C(8) is calculated using the new X in Eq. C(10). 


4) Mach number is calculated.using Eq. B(4). 
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Tale Ci 


Flow Rates Used for Inlet Velocity Determination 
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FLOW FUNCTION 
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APPENDIX D: ROTOR LOSSES 


DI. NASA LOSS COEFFICIENT 

Losses incurred in flows in turbomachinery may be de- 
Scribed with a varity of methods depending upon the measure- 
ments available and the specific application. Ref. 1 


describes a loss coefficient defined by 


fae, - ce 
ws D(1) 
i mt, 
wnere 


/ 

Wo = mass-averaged defect in 
relative total pressure 
divided by the pressure 
equivalent of the inlet 
velocity head 


Ak 
i] 
P, = relative total pressure 
upstream 
a 
>, = relative total pressure 
downstream 





Py = static pressure upstream 


¢ 
(Po)s4 = ideal rebhative total 
pressure downstream 


xj 


rom the measurements obtained and a knowledge of the 


inlet velocity profile the parameters Pi Pos Fe Plos Xyy 


Xx Lee: pe 


—1 
a t2 are known. A representation of W in terms 
of measured stagnation conditions and velocities (X or M) 


is desired, 
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Di,i Upstream of the Rotor, 


may be written: 





The following relations 


T 
= me = 2 ym ) 
Faia ~ Pty Sy? = Ey [t= XG | anes 
we Us - US 
‘se MY ac, 205 D(3) 
T Z 
7 Z 2 2 
TL, ee Kay + Xo - XU D(4) 


where X is a velocity non-dimensionalized as in Eq. B(1), 


Subscript 


denotes the wheel Speed, 


denotes the relative velocity, and subscript W 


The other notation is given in 


the T-S diagram for the flow through the rotor shown in the 


above sketch, 
Using the above 
Por whieh theme is a 


relations result: 


Ey 
ao 1+ XE 
ol 2 
Py = Poy ( 1 + 


and 


Di1.2 Downstream 


exist at the machine 


Kor Xo, and Up where 


equations, for an “impulse” machine 


meridional entrance flow, the following 


D(5) 

SP COR es2Y) 
Xo, ) D(6) 
Ke ) Wy <1 ) D(7) 


of the Rotor, Sufficient measurements 
station 2 to derive 4. using Tt, Tos 


ee is the measured downstream velocity, 


non-dimensionalized using the downstream limiting velocity 


Lae 





(Veo) in Eq. B(i). Alternatively, since T,= constant for 
the flow through the rotor (see diagram), then one of the 
Separate measurements is not needed, Two methods were 
therefore used. Method 1 used all measurements at station 
2. Method 2 did not use the temperature rise at station 2, 
and thus served as a cneck on the accuracy of the temper- 
ature measurement, 

(i) Method I (all measurements at station 2) 


T ¥ 
Pas 2 RD) A8-1) D(8) 


o 
where a is given b 
T+. g ‘i 


pb T+ = 4 [T+ 
BEo=1. xé c=) = 2X, Xu. :. sin, D(9) 


Tto U2 T+2 
Uetne Ba. DCG), D(7), D(8), and D(9) in Eq. D(1), there 





results 


” 


Vey 2XXy upsireela? me 
a [ 1+XG,] PON ore i aa T+o/T +1 Fa) 
é E spe KE alia 2 2 2 xe | ¥/($-1) 


U 
; D(10) 


-(ii) Method II (Temperature measurement at station 2 
not needed) 


Temperature measurement downstream may be found using 


tne following relation. 


TS 
Ev Tr /(Tte/TH1) = (14X 9 )/(Teo/Tey) Dirt) 
Define T+5 ; 2 
-_* X oXUp sind, +\/(X:Xyp sin ,)" + 4 D(12) 


pat 





Semreesa. DVS), D(7), D111), and ™D(12) in Ea. D(1), the 


loss coefficient is given by 


27 (4-1) 2 21 Yy-1) Pte 
| vd [ ‘ U2)/Y Se Cs) 





fe x2 V4) [i - x2) Y(¥-1) 


Peay, LDMAPUSTON PACTOR 

Ref, 1 reports results of correlating compressor losses 
my terms Of the *diffuSion factor". This approach results 
from the argument that the loss in stagnation pressure de- 
pends largely on the suction-Ssurface adverse pressure gra- 
dient, of which tne diffusion factor can be Shown to be an 
approximate measure. The diffusion factor of Ref. 1 is 


defined by 
W, RWy4 ~ RoWys 


W 


SS — D(14) 
fl Oy Sane), 


D=1- 


In thas fomm, changes in ax@#al velocity and radius through 
the blade row are accounted for, 

It is concluded in Ref. 1 that the total-pressure-loss 
parameter ("2982 2) can be correlated in terms of the 
above diffusion factor. Each family of blade shapes has a 
somewhat different correlation curve, however, this pre- 
sents no real problem as the total number of shapes is 


limited, The correlation for several types of blading is 


shown dim bie, 203 %0f Ref. 1. 
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